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ABSTRACT 


In  this  thesis,  an  investigation  was  performed  to 
analyze  the  dynamic  stability  characteristic  of  an  aircraft 
which  has  sustained  damage  to  a  primary  control  surface.  The 
analysis  was  performed  using  the  existing  functional  form  of 
actual  wind  tunnel  data  taken  on  a  F-16  model.  Two  control 

i 

schemes  are  used  for  trimming  an  F-16  that  has  sustained 
damaged  to  its  rudder.  The  First  control  scheme  represent 
the  basic  aircraft,  while  the  second  allowed  the  Horizontal 
Tail  Ailerons  to  move  independently  from  the  Flaperons. 

The  investigation  was  conducted  for  one  flight 
condition  representative  of  the  aircraft  at  cruise  speed. 
Region  in  a//3  space  where  trim  can  be  achieved  was  selected 
as  input  into  a  linearized  aircraft  model.  This  model  took 
into  account  the  failed  control  surface.  The  eigenvalues  of 

the  open  and  closed  loop  models  were  analyzed  to  determine 

) 

the  region  in  a//3  space  where  the  aircraft  was  dynamically 

« 

stable.  The  migration  of  the  eigenvalues  for  several  trim 
conditions  was  also  investigated  to  gain  some  insight  on  the 
aircraft  behavior  while  in  an  unsymme trical  orientation. 

For  this  study,  che  open  loop  eigenvalues  for  the  trim 
area  investigated  gave  a  stable  system.  When  the  aircraft 
controller  was  added  into  the  system,  regions  of  dynamic 
instability  appeared.  For  Rudder  Failure  less  than  -20 


degrees,  trim  could  be  achieved  but  the  aircraft  was 
dynamically  unstable. 
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DYRAMIC  A1ALTSIS  OF  COMBAT  AIRCRAFT 
KITH  COITROL  SURFACE  FAILURE 

I.  IVTRODUCTIOM 


In  modern  high  performs  aircraft.  Flight  Control 
Systems  (FCS)  are  critical  in  achieving  the  performance 
levels  and  operational  utility  required.  Also,  new  designs 
which  increase  the  performance  make  the  aircraft  more 
dependent  on  the  FCS  for  stabilization.  If  a  control  surface 
is  damaged  or  not  operational,  the  control  laws  designed  for 
the  healthy  aircraft  cease  to  be  valid  since  any  signal 
going  to  the  damaged  control  surface  will  be  ignored. 

Studies  showed  that  the  FCS  contributed  up  to  20%  of  the 
aircraft  losses  in  combat  [l«lj.  The  principal  reasons  were 
the  physical  damage,  the  loss  of  function,  or  seriously 
degraded  flying  qualities. 

In  recent  years,  several  methods  have  been  examined  to 
address  the  problem  of  damaged  or  failed  control  surfaces. 
The  development  of  techniques,  like  restructuring  the  FCS, 
to  restore  control  may  have  major  implications  in  aircraft 
flight  safety  ,  sortie  generation  in  a  combat  environment, 
in  reliability  and  maintainability,  and  in  saving  the 
pilot's  life  and  the  aircraft.  Before  considering  applying 
any  of  these  techniques,  we  must  understand  the  dynamics  of 


1-1 


the  aircraft  and  evaluate  whether  stabilization  is  possible. 


Problem  Definition 

If  a  control  surface  is  damaged  or  inoperable,  several 
negative  effects  will  be  encountered.  First,  any  input  going 
to  the  damaged  control  surface  will  be  ignored.  The  FCS  will 
have  to  rely  on  the  other  surfaces  to  control  the  aircraft 
attitude.  Second,  the  coupling  effect  between  the 
longitudinal  and  lateral  modes  of  the  aircraft  may  not  be 
r*.gligible .  For  example,  if  the  rudder  fails  and  is  locked 
into  a  position  other  then  zero,  the  aircraft  is  likely  to 
experience  unwanted  lateral  force  as  well  as  yaw. 

The  questions  that  arise  when  a  control  surface  becomes 
inoperable  are:  could  we  maintain  the  aircraft  in  an 
equilibrium  or  trimmed  state,  and  is  the  aircraft 
dynamically  stable?  Depending  on  the  flight  condition,  many 
newer  types  of  combat  aircraft  have  to  rely  on  the  FCS  to 
provide  dynamic  stability  even  in  a  trimmed  state. 

This  research  will  deal  with  the  latter  question  and 
will  attempt  to  provide  a  better  understanding  of  the 
problem  and  the  means  available  to  address  it. 


Previous  Work 

Eslinger  [2]  investigated  a  failure  of  the  AFTI/F-16 
right  horizontal  tail  with  all  other  surfaces  operational. 
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The  failed  surface  was  left  free  floating.  His  model 
utilized  constant  aerodynamic  derivatives  at  the  selected 
flight  conditions.  In  order  to  restructure  the  control  laws 
for  both  the  healthy  and  damaged  aircraft,  he  used  the 
multivariable  design  technique  developed  by  Professor 
Porter.  As  he  noted  [1],  the  left  horizontal  tail  assumes 
primary  pitch  control  while  the  other  surfaces  deflect  to 
counter  the  rolling  and  yawing  moments  produced  by  the  left 
horizontal  tail  deflection.  Weiss  et  al,  [3],  developed  and 
solved  an  automatic  trim  problem  for  restructurable  aircraft 
control.  In  their  paper,  the  failure  is  treated  as  a 
disturbance  from  desired  steady-state  outputs.  Using  the 
observable  part  of  those  disturbances  that  exist  after  a 
control  surface  failure,  they  feed  forward  a  control 
solution  which  is  a  function  of  the  desired  steady-state 
output  and  the  observed  disturbance.  They  also  noted 
[3:405],  that  the  most  challenging  single  element  failure  is 
a  stuck  rudder  since  it  is  used  extensively  for  damping  the 
dutch  roll  mode,  and  little  side  force  can  be  produced  by 
the  other  control  surfaces. 

Thural,  [4],  conducted  wind  tunnel  experiment  to 
investigate  the  effect  of  various  types  of  control  surface 
failures  on  the  aircraft  stability  derivatives.  He  conducted 
his  test  on  a  one-twentieth  scale  model  F-16  in  the  AFIT 
five  foot  wind  tunnel.  He  collected  data  for  three  different 
configurations,  where  each  represented  a  potential  failure 
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type.  The  data  was  collected  by  varying  each  control  surface 
individually  for  a  given  angle  of  attack  (a)  and  sideslip 

angle  (/S).  Therefore,  the  data  includes  information  about 

the  coupling  of  the  static  aerodynamic  stability 

derivatives.  His  experimental  setup  did  not  permit  him  to 

collect  data  for  dynamic  stability  derivatives. 

In  1989,  Zaiser,  [5],  reduced  the  data  collected  by 
Thural  for  one  particular  failure.  He  employed  a  least 
square  curve  fitting  technique  to  develop  polynomial 
functions  which  describe  the  aircraft  static  stability 
derivatives.  After  deriving  the  equilibrium  equations  for 
rectilinear  flight  in  terms  of  the  static  stability 
derivatives,  he  analyzed  the  impact  of  an  actuator  failure 
of  the  rudder  for  the  F-16  aircraft.  He  also  investigated 
different  control  implementations  which  allowed  for  greater 
independence  of  movement  among  the  undamaged  control 
surfaces.  Region  in  the  a/p  space  where  equilibrium  was 
achievable  were  investigated  at  two  different  flight 
conditions . 

At  the  conclusion  of  his  thesis,  Zaiser  made  several 
recommendations  for  follow-on  work,  [ 5  *  6 6  J .  He  stated  that  a 
dynamic  analysis  should  be  performed  to  evaluate  how  the 
aircraft  would  respond  if  trimmed  in  an  unsymmetrical 
orientation . 
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Purpose 


This  research  will  investigate  the  dynamic  stability 
characteristic  of  an  F-16  aircraft  that  has  sustained  damage 
to  its  rudder  actuator.  Static  aerodynamic  coupling  that 
results  from  unsymmetrical  trim  orientation  will  be  included 
in  deriving  a  linear  state-space  model  of  the  aircraft.  For 
a  given  flight  condition,  several  trimmed  conditions  will  be 
investigated  for  dynamic  stability.  The  impact  of  the 
failure  will  also  be  investigated  for  the  current  aircraft 
control  laws. 

Approach 

To  accomplish  the  stated  purposes  of  this  research, 
specific  tasks  are  accomplished  and  presented  in  the 
different  sections  of  this  thesis.  The  force  and  moment 
coefficients  that  Capt  Zaiser  reduced  into  functional  form 
are  used  in  conjunction  with  the  equilibrium  equations  to 
find  a  trim  condition  for  a  specific  control  implementation. 
The  static  stability  derivatives  are  linearized  for  each 
static  equilibrium  condition  and  included  in  the  aircraft 
plant  model. 

The  linearized  equations  of  motion  are  derived  and 
analyzed  to  relate  the  impact  that  the  different  stability 
derivatives  have  on  the  model.  As  Weiss  pointed  out, 

[5:405],  the  actuator  failure  of  the  rudder  is  assumed  to  be 
the  most  significant  single  primary  control  failure.  This  is 
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taken  into  account,  and  the  actual  control  laws  of  the  F-16 
aircraft  are  modified. 

Presentation 

The  analysis  performed  in  this  thesis  is  presented  in 
the  following  chapters.  Chapter  II  gives  an  overview  of  how 
the  equilibrium  regions  are  obtained  for  the  different 
control  implementations  assumed.  The  derivation  of  the 
linearized  equations  of  motion  and  the  formulation  of  the 
F-16  plant  will  be  included  in  Chapter  III.  Chapter  IV  will 
look  at  the  dynamic  stability  of  the  aircraft  for  the 
different  trimmed  conditions.  The  results  of  this  analysis 
are  presented  and  discussed  in  Chapter  V  and  Chapter  VI 
contains  a  summary  of  the  results  of  this  research  and 
recommendations  for  further  study. 
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II.  TtIH  DETERMIIATIOM 


Introduction 

The  analysis  performed  In  this  thesis  is  based  on  data 
obtained  by  Zaiser,  [5],  in  his  Master's  thesis  in  1989.  In 
this  chapter,  a  short  description  of  the  F-16  is  given  along 
with  a  discussion  of  the  results  obtained  by  Zaiser.  More 
specific  trim  conditions  are  also  evaluated  using  techniques 
similar  to  those  employed  by  Zaiser. 

Aircraft  Description 

The  F-16  is  a  single  engine,  low  aspect  ratio  fighter 
aircraft  currently  in  service  with  several  countries.  Seven 
control  surfaces  are  employed  on  the  aircraft.  All  seven 
control  surfaces  are  of  interest  in  this  research.  The 
location  of  each  control  surface  can  be  found  by  referring 
to  Appendix  A. 

The  primary  function  of  the  Leading  Edge  Flaps  (LEFs) 
is  to  vary  the  camber  of  the  wing  as  the  angle  of  attack  (a) 
increases.  This  causes  CTi|ni#  to  occur  at  higher  a,  thus 
providing  more  lift.  They  are  designed  to  deflect 
symmetrically  and  their  deflection  is  scheduled  as  a 
function  of  a  and  Mach  number.  Therefore,  the  pilot  has  no 
direct  control  authority  on  their  deflection. 

The  Flaperons  (FLs)  are  used  to  provide  both  lift  and 
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rolling  moment.  Below  a  specific  dynamic  pressure  (q),  the 
FLs  act  as  flaps  to  provide  lift.  Otherwise,  they  act  as 
ailerons  which  are  controlled  by  the  pilot  to  provide 
rolling  moment. 

The  Horizontal  Tails  ( HTs )  are  employed  as  elevators  to 
provide  a  pitching  moment  commanded  by  the  pilot.  The  HTs 
also  deflect  asymmetrically  to  augment  the  rolling  moment, 
and  their  deflection  is  scheduled  as  a  function  of  altitude. 


q,  and  Flaperons  input. 

The  Rudder  is  the  primary  control  surface  for  yawing 
the  aircraft.  On  the  F-16,  the  rudder  is  the  dominant 


surface  for  generating  side  forces. 


Table  2-1  F-16  Reference  Data 


Gross  Weight  gw 

Wing  Area  S 

Span  b 

MAC  C 

Center  of  Gravity  Cg 


21018  lbf 
300  Ft2 
29  Ft 
10.94  Ft 
0.35  MAC 


Moment  of  inertia  in  Body  axis 

X  Moment  1^  10033.43  Slug  Ft2 

Y  Moment  Iyy  53876.27  Slug  Ft2 

Z  Moment  Izz  61278.45  Slug  Ft2 

X-Z  Moment  IX2  282.13  Slug  Ft2 

Control  Surface  Deflection  Limits 


LEF 

-2° 

a 

6 

a 

25 

FL 

-20° 

s 

6 

a 

20 

HT 

-25° 

a 

6 

a 

25 

Rudder 

-30° 

a 

6 

a 

30 

2-2 


The  sign  convention  adopted  in  this  thesis  is  shown  in 
Figure  A-2  of  Appendix  A.  Basic  aircraft  data  is  presented 
in  Thural  thesis  [4«27],  and  is  summarized  in  Table  2-1. 


Aerodynamic  Coefficients 

The  data  collected  by  Thural,  [4],  is  presented  as 
nondimensional  force  and  moment  coefficients  for  a  given  a 
and  j9 .  Each  coefficient  can  be  transformed  into  a  force  or 
moment  using 


Pg  -  CF  Q  S 
Nm*  CHqS  b 


(2-1) 


where  Fg  and  Ms  represent  the  appropriate  force  or  moment 
acting  on  the  aircraft  in  the  stability  axis  system.  Figure 
A-2  gives  a  graphical  representation  of  each  axis  system. 

For  a  rigorous  definition,  refer  to  Etkin  work  [7:106-112]. 
Since  the  data  was  taken  at  finite  discrete  points,  it  was 
transformed  in  functional  form  for  analytical  purpose.  Using 
a  Least  Square  curve  fitting  technique,  Zaiser  determined 
the  contribution  of  each  force  and  moment  on  the  aircraft 
and  its  control  surfaces  as  a  function  of  a  and  p,  [5:10- 
16].  The  functional  form  of  each  force  or  moment  coefficient 
are  represented  by  equation  2-2  where  the  first  term 
represents  the  contribution  of  the  basic  aircraft  with  no 
control  surface  deflections,  and  the  second  term  represents 
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Equilibrium  State 


The  desired  trim  conditions  that  are  investigated  in 
this  research  have  the  aircraft  flying  a  rectilinear 
trajectory  at  constant  altitude.  When  the  aircraft  is  in 
equilibrium,  all  external  forces  and  moments  acting  on  it 
equal  zero.  Therefore  the  equilibrium  equations  become 


+  FTx  -  mg  sinfl  =  0 

(2-3) 

+  mg  cos©  sin$  =  0 

(2-4) 

*Vs 

+  mg  cos0  sin$  =  0 

(2-5) 

M»  =  0 

(2-6) 

•< 

H 

O 

(2-7) 

II 

O 

(2-8) 

where  FAi,  F^,  and  represent  the  aerodynamic  forces,  the 
thrust  force,  and  the  aerodynamic  moments  in  the  i  axis  of 
the  body  axis  reference  system.  9  and  ♦  are  the  Euler  angles 
that  define  the  aircraft  attitude  with  respect  to  earth 
inertial  reference  frame.  The  only  other  equation  required 
in  the  trim  analysis  is  an  expression  that  defines  the 
aircraft  pitch  angle  for  a  constant  altitude  flight,  which 
is 

0  =  Tan'1  |  Tana  Cos$  +  cewi  )  <2_9> 
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Since  no  restriction  was  placed  on  a  wing  level  flight, 
coupling  effects  between  the  longitudinal  and  lateral  Bodes 
are  apparent  in  equations  2-4,  2-5  and  2-9.  A  complete 
derivation  of  the  equilibrium  equations  can  be  found  in 
[5.85-99] . 

Before  solving  the  trim  equations  for  a  specific  rudder 
failure,  flight  conditions  need  to  be  established.  Table  2-2 
gives  the  flight  condition  that  is  investigated  in  this 
research . 


Table  2-2  Flight  Condition 


Mach 

0.6 

Altitude 

15000  Ft 

Velocity 

375  KEAS 

q 

300  psf 

Problem  Formulation 

The  desired  equilibrium  is  a  rectilinear  flight  at 
constant  altitude.  Although  equilibrium  states  might  be  less 
difficult  to  achieve  at  other  flight  conditions,  only 
rectilinear  flight  will  be  investigated.  A  failure  of  the 
rudder,  which  results  in  the  rudder  being  locked  into  a 
specific  deflection  is  the  only  failure  mode  that  this 
thesis  will  study.  The  investigation  will  also  be  limited  by 
the  range  of  the  test  data  that  was  collected  by  Thural  [4]. 
Therefore  the  dimensions  of  the  a/p  space  that  will  be 
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investigated  are  limited  to  -6.0°  s  p  a  6.0°  and  0°  s  a  s 

20.0°. 

Some  assumptions  still  need  to  be  stated  before 
proceeding  with  the  analysis.  They  aret 

1.  The  aircraft  is  assumed  to  be  a  rigid  frame. 

2.  The  earth  surface  is  assumed  to  be  an  inertial  frame 

of  reference. 

3.  The  aircraft  mass  and  mass  distribution  are  assumed 

to  be  constant. 

4.  The  X-Z  plane  of  the  aircraft  is  assumed  to  be  a 

plane  of  symmetry. 

These  assumptions  hold  for  both  the  equilibrium  equations 
and  the  linearized  equations  of  motion  developed  in 
Appendix  D. 

As  discussed  in  the  beginning  of  this  chapter,  the  only 
control  authority  that  the  pilot  has  on  a  healthy  aircraft 
is  through  the  Horizontal  Tail  Elevator  (HTE),  the  Rudder, 
and  the  Flaperons  (FLs).  The  Horizontal  Tail  Ailerons  (HTAs) 
deflection  is  proportional  to  the  FLs  deflection.  For  the 
flight  conditions  of  Table  2-2,  the  HTAs  deflect  only  a 
factor  of  0.294  of  the  FLs  {8]. 

The  control  schemes  investigated  in  this  thesis  are 
derived  by  allowing  successively  greater  independence.  One 
point  to  note  is  that  the  control  scheme  discussed  in  this 
chapter  does  not  refer  to  the  control  laws.  The  two  control 
schemes  investigated  are  shown  in  Table  2-3. 


2-7 


Table  2-3  Control  schemes 


Case  A 

Case  B 

6FL 

6fl 

6hte 

6hte 

6  HTA 

Case  A  represents  the  basic  aircraft.  Case  B  allows  the 
HTA  to  deflect  independently  from  the  FLs.  It  is  assumed 
that  an  algorithm  is  present  to  trim  the  aircraft  using 
these  control  surfaces  without  modifying  the  actual  control 
laws.  Also,  the  deflections  of  the  individual  control 
surfaces  in  both  cases  are  related  as  follows! 

&fl  ~  ^  (&xn  ~  &un)  (2-10) 

^BTE  =  ~2  ^XBT  +  ^  lot)  (2~11) 

^  HTA  ~  ~2  ^ XBT  ~  ^  LOT1  (2~12) 

=  \  +  <2'13) 


Before  solving  the  trim  problem,  the  external  forces 
and  moments  acting  on  the  aircraft  still  need  to  be 
determined.  By  specifying  a,  0,  and  the  dynamic  pressure  q. 
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the  external  forces  and  moments  are  evaluated  using  the 
polynomial  listed  in  appendix  B  as  a  function  of  a  specific 
control  surface  deflection.  The  total  external  force  or 


mon.jnt  acting  on  the  aircraft  can  be  written  as 


Ft  -  A.  ♦  s  *£  £  £  c,«  «”  P* 

1^1  m-0  n-0 


(2-14) 


where  F^  is  the  total  force  or  moment  acting  on  the 
aircraft,  A0  is  the  contribution  of  the  aircraft  with  all 
the  control  set  to  zero,  B  is  the  contribution  of  the  failed 
rudder  and  the  LEFs,  and  the  last  term  is  the  force  or 
moment  that  results  from  the  unknown  deflection  of  the 
control  surfaces.  The  unknowns  th*»t  regain  to  be  evaluated 
are  the  deflection  of  the  control  surface. 

Solving  the  trim  problem 

Assuming  that  the  power  available  from  the 
aircraft  can  compensate  for  the  aerodynamic  forces  and  the 
gravitational  term  in  equation  2-3,  the  problem  can  be 
formulated  as  follows 

i 

-  ( Az  +  Bz  +  mg  CosO  Sin+)  =  g  (2-15) 

-  (Au  +  Bm)  =  g  (2-16) 
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(2-17) 


-  (al  +  bl)  =  g  6, 

+  Bh)  =  ^ 


Since  the  components  on  the  left  side  of  each  equation  are 
known,  the  problem  can  be  reformulated  as 

b  -  [  A  ]  6  (2-19) 

where  the  vector  b  represents  the  known  forces  and  moments, 
the  A  matrix  contains  the  control  derivatives,  and  6  is  the 
unknown  control  deflection  vector.  Using  equations  2-4,  2-9, 
and  2-15  to  2-18,  the  trim  problem  can  be  solved.  For  Case 
A,  the  I  in  equations  2-15  through  2-18  is  two,  since  only 
the  FLs  and  the  HTE  are  directly  controlled  by  the  pilot.  In 
case  B,  I  equals  three  since  the  HTAs  are  assumed  to  be 
independent  from  the  FLs. 

Equation  2-4  and  2-9  are  used  first  to  estimate  4  and 
8.  For  both  cases,  equation  2-19  represents  an 
overdetermined  system  of  equations  that  can  be  solved  using 
Singular  Value  Decompositions  (SVD)  [9*59]. 

In  his  thesis,  Zaiser  wrote  a  computer  code  to  solve 
for  the  trim  conditions,  where  the  A  matrix  in  equation  2-19 
was  square  [ 5 1 1 1 2 ] .  The  same  code  is  modified  and  used  to 
evaluate  the  trim  conditions  for  each  case  at  a  given  rudder 
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failure.  Zaiser  program  was  modified  by  the  inclusion  of  a 
least  square  algorithm  in  the  SVD  routine.  Since  a  least 
square  solution  can  be  obtained  for  any  given  set  of  initial 
conditions,  a  range  of  values  that  would  be  considered  zero 
had  to  be  established.  Table  2-4  gives  the  limits  that  are 
incorporated  into  the  code. 


Table  2-4  Forces  and  Moments 

Limits 


1  FAzi 

< 

50.0 

lbf 

1 

< 

50.0 

lbf 

1  MAxl 

< 

500.0 

lbf 

Ft 

lMAyl 

< 

500.0 

lbf 

Ft 

1  mAz^ 

< 

500.0 

lbf 

Ft 

Results  from  the  Trim  Analysis 

The  results  from  the  trim  analysis  are  presented  in 
Figures  2-2  to  2-7.  The  area  where  equilibrium  is  possible 
is  presented  on  the  first  graphic  of  each  figure  for  a  given 
case  and  rudder  failure.  The  roll  angle  and  the  control 
surface  deflections  where  the  aircraft  is  in  equilibrium  are 
also  presented.  For  each  of  the  cases  and  rudder  failures, 
if  /3  is  specified,  only  a  very  small  variation  in  a  is 
permitted  for  the  aircraft  to  remain  in  equilibrium. 


2-11 


TRIM  AREA  „  *  ROLL  ANGLE 


(93G)  330NV  TOU 


rn  rd  ro  r*S 


o 

W 

Q 

£ 

w 

CQ 


(030)  VHdlV 


£ 

O 

P 

u 

►J 

bn 

w 

Q 

hJ 

e 


(03a)  NOII03333Q  TttLO 


Figure  2-2  Trim  Region  for  Case  A 
0  Degree  Rudder  Failure 
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Figure  2-3  Trim  Area  for  Case  A 
-10  Degrees  Rudder  Failure 
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Figure  2-4  Trim  Area  for  case  B 
0  Degrees  Rudder  Failure 
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-10  Degrees  Rudder  Failure 
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Figure  2-6  Trim  Area  for  Case  B 
-20  Degrees  Rudder  Failure 
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Figure  2-7  Trim  Area  for  Case  B 
-25  Degrees  Rudder  Failure 
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For  a  zero  degree  rudder  failure,  both  cases  show  some 
asymmetry.  This  results  from  the  limits  of  Table  2-4  that 
needed  to  be  included  in  the  code  to  calculate  the  trim 
area.  For  a  rudder  failure  other  than  zero  degrees,  a 
vertical  line  is  included  in  each  graphic  at  0«-6°  since  the 
functional  representation  of  the  stability  derivative  is 
only  valid  for  -6°  <  £  <  6°  and  0°  <  a  <  20°,  [5i30].  Data 
points  to  the  left  of  that  line  are  not  analyzed,  since  for 
p  <  -6°,  the  results  are  obtained  by  extrapolating  the 
functional  form  of  the  stability  derivatives. 

One  of  the  interesting  features  displayed  in  the 
results  is  the  attitude  of  the  aircraft  for  a  specific 
condition.  For  a  25°  rudder  failure,  case  B,  it  might  be 
preferable  to  trim  the  aircraft  at  j8*-6°,  a-1.8°  and  *--18° 
since  this  is  the  trim  condition  that  gives  the  most  control 
authority  to  each  control  surface.  This  will  be  taken  into 
consideration  in  the  following  chapters. 

Summary 

In  this  chapter,  different  rudder  failure  are  presented 
for  a  given  control  scheme  and  flight  condition.  The 
aircraft  controls  are  also  presented  with  their  limitations. 
Assumptions  are  made  regarding  the  control  scheme 
implementation  for  both  cases.  Data  required  to  analyze  the 
aircraft  dynamic  response  for  a  specific  failure  conditions 
is  presented. 
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III.  AIRCRAFT  PLANT  MODEL 


Introduction 

In  this  chapter,  the  linearized  F-16  plant  model  is 
derived  using  the  functional  representation  of  the  static 
stability  derivatives  of  Appendix  B  for  several  trim 
conditions.  The  eigenvalues  of  the  open  loop  plant  are 
analyzed  to  determine  the  dynamic  characteristics  of  the 
aircraft  for  a  specific  rudder  failure. 

Linearized  equations  of  motion 

The  mathe’-  -  cal  model  is  presented  in  Appendix  D.  The 
sign  conven.i  ,1  for  the  axes  and  control  deflections  are 
shown  in  Figure  A-2  of  Appendix  A.  The  simplifying 
assumptions  made  during  the  derivation  of  the  equations  of 
motion  are: 

Assumption  1.  The  aircraft  is  assumed  to  be  a  rigid 

body . 

Assumption  2.  The  earth  surface  is  assumed  to  be  an 

inertial  reference  frame. 

Assumption  3.  The  mass  and  mass  distribution  of  the 

aircraft  is  assumed  constant. 

Assumption  4.  Disturbances  from  steady  flight 

conditions  are  small,  implying  a  small 
angle  approximation.  Higher  order  terms 
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of  the  disturbance  quantities  are 
negligible . 

Assumption  5.  The  flow  surrounding  the  aircraft  is 

assumed  quasi-steady. 

Assumption  6.  Variation  of  the  atmosphere,  including 

density  and  speed  of  sound,  are 
negligible  for  small  altitude 
perturbations . 

Even  if  these  assumptions  are  made,  the  equations  still 
include  coupling  effects.  The  complexity  can  be  reduced  by 
linearizing  the  equations  of  motion  about  a  steady  state 
flight.  For  this  investigation,  the  only  restriction  is  for 
the  aircraft  to  fly  a  rectilinear  trajectory  at  constant 
altitude  which  implies! 

1.  Initial  side  velocity  may  exist:  VQ 

2.  Initial  bank  angle  may  exist  : 

3.  Initial  pitch  angle  exists  :  eQ 

4.  Ho  initial  angular  velocities  exist  which  results  in 


P0  =  £0  =  i?0  =  *  =  0  =  4>  =  0  (3'1) 

Equation  3-1  considerably  simplifies  the  linearized 
equations  of  motion  shown  in  Appendix  D,  but  they  still 
include  coupling  effects  produced  by  the  roll  angle  and  the 
side  velocity.  In  this  research,  the  effects  of  unsteady 
atmospheric  disturbances  were  neglected. 
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Aerodynamic  Forces  and  Moments 


Assuming  that  the  aircraft  has  enough  thrust  so  that 
steady  state  flight  can  be  maintained,  the  thrust  setting 
will  balance  the  remaining  forces  in  the  corresponding  axis. 
Also,  the  direct  thrust  contributions  to  the  stability 
derivatives  is  generally  negligible  for  conventional 
aircraft  and  is  assumed  to  be  zero  for  this  aircraft 
[6:267].  The  only  part  of  the  linearized  equations  of  motion 
that  still  need  development  are  the  aerodynamic  forces  and 
moments  . 

The  representation  of  the  aerodynamic  forces  and 
moments  is  usually  made  in  the  stability  axis  system.  Since 
the  perturbed  equations  of  motion  are  written  in  the  body 
axis  system,  the  aerodynamic  forces  and  moments  need  to  be 
transformed  into  that  axis  system.  Each  stability  axis  can 
be  transformed  into  the  body  axis  system  using  the  stability 
to  body  axis  transformation  matrix  ( [ BS  ]  )  presented  in 
Appendix  D.  Once  the  transformation  is  made,  the  aerodynamic 
forces  and  moments  are  given  in  the  same  axis  system  as  the 
equations  of  motion. 

Each  of  the  forces  and  moments  need  to  be  expanded  to 
determine  their  dependence  on  the  perturbed  motion.  The 
expansion  is  done  using  a  Taylor  series  expansion  at  a  given 
trim  condition,  denoted  by  the  subscript  ()0.  The  expansion 
of  the  forces  and  moments  can  be  represented  as 
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F‘  F“  *  (m;)1"*1  *  *•••*  (l^)l0A» 

M  =  w°  +  (J^)I°X1  ‘•••t('^)l°A” 


where  the  variable  A  represent  the  perturbed  velocities  and 
accelerations,  and  X  the  perturbation  from  the  ()0 
condition.  The  higher  order  terms  have  been  eliminated  from 
the  expansion  in  accordance  with  assumption  4.  The  results 
of  this  expansion  is  shown  in  Appendix  D.  The  analysis  of 
the  aircraft  motion  is  performed  at  different  trim 
conditions  where  the  aircraft  attitude  preclude  the 
separation  of  the  longitudinal  and  lateral  mode.  Therefore, 
all  cross-coupling  derivatives  are  included  in  the  perturbed 
equations  of  motion. 

State  Space  Form 

In  order  to  analyze  the  system,  the  perturbed  equations 
of  motion  are  put  into  matrix  form.  This  form  is  used  to 
determined  the  dynamic  stability  and  control  for  various 
rudder  failure  and  trim  conditions.  The  matrix  form  of  the 
perturbed  equations  of  motion  may  be  written  asi 

X  =  A  x  +  B  u  (3~4> 
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where  x  represent  the  aircraft  states,  A  the  plant  matrix,  u 
the  control  variables,  and  B  the  matrix  of  coefficients 
associated  with  the  control.  To  determine  the  stability  of 
the  system,  the  eigenvalues  of  the  plant  matrix  are  computed 
using  a  control  analysis  program  (Pro-Matlab  [10]).  If  an 
eigenvalue  has  a  positive  real  part,  then  this  state  is 
unstable.  Dyn.  uic  stability  still  might  be  recovered  with 
the  controller.  This  aspect  will  be  discussed  in  the  next 
chapter.  Each  of  the  eigenvalues  are  also  presented 
graphically  to  gain  some  insight  on  the  plant  behavior.  It 
is  also  possible  to  analyze  the  controllability  of  the 
system.  Equation  3-5  is  one  method  to  determine  the 
controllability  of  the  system.  In  this  equation,  A  is  the 

Mc  =  [  B  5  AB  :  A*B  i  ...  I  Aa~xB  ]  (3-5) 


plant  matrix,  B  the  control  matrix,  Mc  the  controllability 
matrix,  and  n  the  number  of  states.  If  Mc  is  full  rank, 
rank(Mc)  =  n,  then  that  particular  trim  condition  is 
completely  controllable  with  the  available  inputs  using 
state  feedback.  If  the  system  is  completely  controllable 
then  it  is  possible  to  reach  any  state  [11:2-42]. 

Results 

In  order  to  facilitate  the  data  handling,  a  computer 
code  was  written  to  determine  the  plant  matrix.  A,  and  the 
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control  matrix,  B,  for  different  trim  conditions.  The 
computer  code  is  presented  in  Appendix  E.  The  static 
stability  derivatives  of  Appendix  B  are  linearized  and 
included  into  the  code.  For  specific  a  and  p  the  static 
stability  derivatives  are  evaluated  and  included  in  the 
linearized  equations  of  motion.  Table  3-1  lists  the  static 


Table  3-1  Static  stability 
Derivatives 


stability  derivatives  that  are  linearized.  The  remaining 
stability  derivatives  were  taken  from  data  collected  by  the 
Flight  Dynamics  Laboratory  for  the  AFTIF-16  flying  at  Mach 
0.6  for  three  different  altitudes  (0,  5000,  and  30000  FT). 
Since  the  altitude  of  interest  in  this  research  is  15000 
feet,  the  stability  derivatives  were  estimated  using  a 
second  order  polynomial  fit  of  the  three  data  point 
available.  Table  3-2  list  the  value  of  the  derivatives  that 
are  included  in  the  computer  code. 

Since  no  data  on  the  dynamic  cross-coupling  derivatives 
were  available  for  the  F-16,  they  were  assumed  to  be  zero. 
This  is  a  reasonable  assumption,  since  for  low  or  their 
contribution  to  aircraft  motion  is  relatively  small 
according  to  Orlik-Rukemann,  [ 1 2  s 1 - 1  ]  .  For  each  of  the  cases 
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Table  3-2  control  derivatives 


CD 

CL 

Cm 

u 

0.000059 

0 . 00002 

-0.0000638 

a 

0.0 

-0. 99333 

-0.77776 

q 

0.0 

2 . 3989 

-2 . 6761 

CY 

C1 

Cn 

fi 

0.0 

0.0 

0.0 

p 

0.080111 

-0. 23708 

-0.0079264 

r 

0. 53755 

0.025172 

-0.48192 

and  rudder  failures,  the  eigenvalues  are  represented 
graphically  in  Figures  3-1  to  3-8.  The  vertical  line  at 
x=0 . 0  on  most  figures  is  a  reference  line.  The  eigenvalues 


for  case  A  do  not  change  much  over  the  trim  area  for  both 
rudder  failures  (less  than  are  2%)  therefore  an  average 
value  of  each  of  them  is  considered  appropriate.  They  are 
presented  in  Table  3-3. 


Table  3-3  Eigenvalues  for  Case  A 


0  Degree  rudder 

failure  -10 

Degree  Rudder  Failure 

3. 16  s  a  s  3 

.  19 

0.88  s  a 

s 

0.95 

0.05  s  p  s  0 

.  10 

-3.28  s  p 

s 

-3.22 

0.19  s  *  s  0 

.  42 

-7.75  s  * 

s 

-7.48 

-0.0023  ±  0 . 06 76 i 

Phugoid 

-0.0022 

t 

0 . 0682 i 

-0.7275  t  4 . 6 30 1 i 

Short  Period 

-0.6499 

± 

4. 6773i 

-.3624  ±  4 . 0805i 

Dutch  Roll 

-0.4097 

i 

3. 8549i 

-1.4703 

Roll 

-1 . 

5711 

-0.0309 

Spiral 

-0.1 

0199 

The  legend  that  describes  each  figure  is  in  Table 
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Table  3-4  Legend  for  Figures  3-1  to  3-8 


a : 

Dutch  Roll 

b: 

Short  Period 

C  i 

Roll 

d: 

Phugoid  &  Spiral 

RUDDER  FAILURES  DEGREES 

0 

-10  -20 

-25 

*1 

0.0° 

0.0°  -4.4° 

-5.0° 

p2 

-3.3° 

-6.0°  -6.0° 

-6.0° 

The  functional  form  of  the  static  stability  derivatives  was 
derived  from  data,  where  -6.0°  <  P  <  6.0°  and  0.0°  <  a  < 
20.0°.  Therefore  the  curve  fitting  used  for  points  outside 
the  p  limits  may  show  odd  behavior.  This  is  the  reason  why 
the  eigenvalues  have  a  strong  departure  for  p  <  -6.0°  and 
for  rudder  failures  greater  than  -10°.  One  other  interesting 
point  is  the  dynamic  stability  of  the  open  loop  system. 

Since  all  the  real  parts  of  the  eigenvalues  are  negative, 
the  open  loop  plant  is  dynamically  stable.  Equation  3-5  was 
used  to  look  at  the  controllability  of  the  system.  Boundary 
points  as  well  as  some  intermediate  points  were  evaluated 
for  controllability.  All  the  points  that  were  checked 
yielded  a  controllable  system.  This  means  that  with  an 
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Figure  3-1  Case  B  0  Degrees  Failure 
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Figure  3-2  Case  B  0  Degrees  Rudder  Failure 
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Real  Axis 


LOCUS  OF  OPEN  LOOP  EIGENVALUES 


Figure  3-3  Case  B  -10  Degrees  Rudder  Failure 
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REAL  AXIS 


LOCUS  OF  OPEN  LOOP  EIGENVALUES 
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LOCUS  OF  OPEN  LOOP  EIGENVALUES 


Figure  3-7  Case  B  -25  Degrees  Rudder  Failure 
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REAL  AXIS 


appropriate  controller,  it  is  possible  to  position  the 
eigenvalues  of  the  system  to  guarantee  dynamically 
stability.  Now  the  question  that  needs  to  be  answered  is: 
Does  the  present  controller  on  the  aircraft  still  adequate 
considering  the  damaged  rudder?  The  next  chapter  will 
provide  a  means  to  answer  this  question. 

Summary  In  this  chapter,  all  the  equations  needed  to 
analyze  the  dynamic  stability  of  the  open  loop  F-16  aircraft 
having  sustained  a  rudder  failure  were  derived  and  put  into 
matrix  form  to  facilitate  the  analysis.  The  results  showed 
that  the  aircraft  remains  dynamically  stable,  but  with 
lightly  damped  phugoid,  short  period  and  dutch  roll  modes.  A 
rigorous  analysis  of  the  results  is  done  in  Chapter  V. 
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IV.  STATE  SPACE  MODEL  DEVELOPMENT 


Introduction 

In  order  to  investigate  if  the  F-16  control  system 
augments  the  aircraft  dynamic  stability  for  given  trim 
conditions  and  specific  rudder  failures,  a  state  space  model 
of  the  aircraft  flight  control  system  was  created.  The 
aircraft  plant  derived  in  Chapter  III  is  for  different  trim 
conditions,  but  at  a  specific  speed  and  altitude.  To  be 
rigorous,  a  state  space  model  would  have  to  be  developed  for 
each  trim  condition,  since  the  model  is  dependent  on  the 
trim  velocity  of  the  aircraft  and  a  trim.  For  the  range  of 
trim  conditions  analyzed,  the  same  control  system  model  is 
used  throughout  the  analysis. 

Plant  Matrix  Development 

In  order  to  construct  a  state  space  representation  of 
the  F-16  control  system,  a  flight  condition  must  be 
selected.  The  control  law  diagram  presented  in  the  F-16 
Software  Mechanization  Document  [8]  is  linearized  about  the 
flight  condition  presented  in  Table  2-2.  Mo  pilot  input  are 
used,  so  all  paths  associated  with  pilots  inputs  can  be 
ignored.  Since  the  horizontal  tail  is  used  to  command  both 
pitch  and  roll  rates,  an  effective  flaperon  deflection  input 
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was  determined  [8].  The  effective  flaperon  deflection  is 


® Fett  =  +  .2946^ 


(4-1) 


where : 

=  Effective  Flaperon  Deflection  (°) 
=  FlaperonDe flection  (*) 

6^.  =  Horizontal  Tail  Deflection  (#) 


This  effective  flaperon  deflection  was  incorporated  into  the 
computer  code  to  calculate  the  B  matrix  for  both  cases  at 
the  given  trim  conditions.  The  effective  flaperon  deflection 
is  on?  '  used  for  the  roll  rate  commands.  Another 
modification  is  also  made  to  the  control  law  diagram.  The 
load  factor  command  is  change  to  pitch  rate  command.  The 
gain  in  the  command  path  of  the  control  law  diagram  has  to 
be  adjusted  in  order  to  convert  load  factor  to  a  pitch  rate 
command.  This  is  done  using  the  steady  Z  axis  acceleration 
as  shown  in  equation  4-2. 


A 


n 


<JVq 

(  57.3  )  (  32.2  ) 


(4-2) 


where : 

An  =  normal  acceleration  at  pilot  station  ( g ) 
q  =  pitch  rate  (°/s) 

V0  =  steady  state  forward  velocity  (ft/s) 
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Figures  4-1  and  4-2  show  the  final  configuration  of  the 
linearized  control  law  for  both  the  longitudinal  and  lateral 
axis.  Since  the  aircraft  has  a  failed  rudder,  the  lateral 
axis  still  needs  to  be  altered  to  represent  it.  This  is  done 
by  removing  all  feedback  paths  that  are  input  into  the 
rudder.  The  final  configuration  of  the  linearized  control 
law  for  the  lateral  directional  axis  is  shown  in  Figure  4-3. 

The  state  vector  used  to  represent  the  aircraft  is 
shown  in  equation  4-3. 

I-[uaPpgr^6^]  (4-3) 


Since  the  commanded  input  is  pitch  rate  instead  of  load 
factor,  the  outputs  of  the  system  available  for  feedback  are 
a,  q,  and  An.  An  is  in  units  of  g's.  The  expression  for  the 
normal  load  factor  at  the  pilot  station  is 


-  *  -  QV0  -  xa4 


(4-4) 


which  can  be  transformed  using  small  angle  approximation 
into 

a,  -[  -v.  <  *  -  «>  *  *.<*] 

where  X,  is  the  distance  from  the  aircraft  eg  to  the 

d 

accelerometer  located  under  the  pilot's  seat.  For  the  F-16, 


32.2 


(4-5) 
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Figure  4-2  Modified  F-16  Lateral  Control  Systea 
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Xa  is  14.0  feet.  The  next  step  is  to  build  the  measurement 
matrices,  C  and  D,  associated  with  the  plant.  Referring  to 
figure  4-1,  4-3,  and  equation  4-5,  the  measurement  matrices 
become  > 
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where  Anl  and  An2  are  i 


Am 


Ana 


a(2'i)  ♦(so)**5--0 

[000010000] 


(4-7) 


with  i  -  1,  2,  3,  ...,  9,  j-  1,  2,  and  a(2,i),  a(5,i), 

b(2,j),  b ( 5 , j )  representing  the  second  and  fifth  row  of 
elements  of  the  A  and  B  matrix  developed  in  Chapter  III. 
Therefore  the  final  open  loop  system  can  be  represented  by 

i  =  A  X  +  B  U  (4-8a) 
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x  =  C  X  +  D  U 


( 4-8b ) 


In  order  to  prevent  implicit  algebraic  equations  while 
deriving  the  closed  loop  system,  the  state  vector,  X,  can  be 
redefined  by  including  the  control  deflection  into  the  state 
and  the  actuator  model  into  the  input.  This  will  then  lead 
to 

£  =  A*  k  ♦  B'  ( 4-8c ) 

y!  =  d  xf  <4~8c) 

where  the  D  matrix  is  part  of  C'. 


Controller  Development 

The  feedback  and  feedforward  paths  shown  in  Figure  4-1 
and  4-3  can  be  expressed  as  a  matrix  in  the  Laplace  domain 
in  terms  of  the  aircraft  inputs  and  outputs  as 


'UT 


’Faff] 


1 . 056 (s+5) (3S+11.25) 
S(S+ 11 . 25) 

4.2 

(s+io) 

0 

0.353 (3+5) (3g+ll . 25) 
(s+1) (s+11.25) 


(4-9) 


-20.99 (s+5) 
s 
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0 

-0.12 
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The  transf ormation  of  the  matrix  from  the  Laplace 
domain  to  the  time  domain  can  be  done  by  transforming  each 
Laplacian  element  into  a  state  space  phase  variable 
canonical  form  [ 13 « 210-215 ] .  The  transformation  in  this  case 
was  done  using  the  command  called  tfm2ss  of  the  control 
analysis  computer  program  called  Pro-Matlab,  [10].  A  minimum 
realization  was  also  performed  on  the  matrix  to  remove 
unnecessary  states.  The  feedforward  (subscript  E)  and 
feedback  (subscript  K)  in  the  time  domain  are  shown  in 
equation  4-10a  to  4-llb  respectively. 

XB  =  At  XB  +  Bm  (4-10a) 

yB  =  CBxB  +  Dtlcmd  =  iiz  (4-ieb) 

AkXk  +  BkX  (4-lla) 


where 


3  Cjc  ¥■  ~  (4-llb) 


cmd  [  Qcnd  PcMtd  ] 


and 


II  *  +  llj 


( 4-  12b ) 
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Closed  Loop  System  Derivation 


The  F- 16  utilizes  negative  input  and  positive  feedback 
in  its  control  law  diagran,  as  shown  in  Figures  4-1  and  4-2. 
This  is  due  to  the  sign  convention  which  defines  a  positive 
deflection  of  the  effective  flaperon  or  the  horizontal  tail 
as  being  trailing  edge  down.  Since  the  computer  program  used 
to  get  the  closed  loop  state  space  representation  uses 
negative  feedback,  the  sign  of  the  C  and  D  matrix  need  to 
be  changed.  The  block  representation  of  the  total  system  is 
shown  in  Figure  4-4.  Using  Eq  (4-lCa)  through  (4-12b),  it  is 
now  possible  to  derive  the  closed  loop  model  of  the 
aircraft . 

Substituting  equation  4-8d  into  4-lla  and  4-llb  gives 

XK  -  Ak  Xk  +  Bk  C  2d.  (4-13a) 

Ux  -  CK  +  Dk  C  xL  ( 4-13b ) 

Placing  equation  4-12b  into  4-8c  and  substituting  4-13a  and 
4-13b  will  lead  to 

i!  =  (Af  +  B'  DKd)  *  ♦  B'  CKxg  +  B'  Ds  6^  (4-14) 
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Combining  equation  4-10a,  4-14,  4-8d,  and  4-13a  into  matrix 
form  yields 


ab  o  0 

•  • 

X* 

= 

B1  Cg  A'  +  B1  Dgd  B'  CK 

2Sf 

0  Bg  d  Ag 

*g 

+ 


B, 

B'  De 
0 


[®  l«d 


(4-15) 


=  [  o  d  o  ] 


(4-16) 


Using  Equation  4-15  and  4-16  it  is  now  possible  to 
determine  the  characteristics  of  the  closed  loop  system. 
Appendix  F  shows  an  example  of  one  open  loop  and  closed  loop 
systems  for  one  trim  condition  of  case  B. 


Summary 

In  this  chapter  the  closed  loop  system  was  derived 
using  a  modified  controller  that  takes  into  account  the 
failed  rudder.  The  eigenvalues  can  then  be  calculated  for 
each  case  and  failure.  The  open  loop  and  closed  loop 
eigenvalues  can  be  compared  to  determine  the  effectiveness 
of  the  controller.  Analysis  of  the  results  is  performed  in 
the  following  chapter. 
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V.  DISCUSSION  Of  RESULTS 


Introduction 

In  the  previous  chapter,  the  closed  loop  plant  for  the 
F-16  aircraft  was  derived  considering  a  damaged  rudder.  The 
model  is  used  to  determine,  for  each  rudder  failure  and  trim 
condition,  the  closed  loop  eigenvalues  of  the  system.  A 
physical  explanation  is  given  on  the  behavior  of  the  system 
and  its  implication  on  the  aircraft  dynamic  stability. 

Eigenvalues  of  the  Closed  Loop  System 

Using  the  controller  developed  in  the  previous  chapter, 
the  eigenvalues  for  each  rudder  failure  and  trim  condition 
are  calculated.  Since  the  linearization  of  the  equations  of 
motion  presented  in  Appendix  D  includes  cross  coupling 
derivatives  and  asymmetric  trim  condition,  the  results 
obtained  from  the  closed  loop  system  are  dependent  on  the 
angle  of  attack  (a),  the  sideslip  angle  (0),  and  the  bank 
angle  (*).  Table  3-4  in  Chapter  III  also  describe  Figures 
5-1  to  5-8. 

Longitudinal  Motion 

For  case  A,  the  closed  loop  eigenvalues  for  a  rudder 
failure  of  0  and  -10  degrees  are  presented  in  Figures  5-1 
and  5-2.  The  controller  in  this  case  does  not  significantly 
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change  the  location  of  the  eigenvalues  close  to  the 
imaginary  axis. 

If  the  rudder  fails  at  zero  degrees  for  Case  B,  Figures 
5-3  and  5-4  show  the  migration  of  the  eigenvalues  from  p^-  0 
degrees,  4^*  0  degrees,  3.17  degrees  to  02«  “3.3 

degrees,  #2«  "5  degrees,  a2=  .8  degrees.  The  phugoid  complex 
conjugate  eigenvalues  remain  very  close  to  the  origin  (d), 
for  both  the  open  and  closed  loop  system,  since  they  depend 
primarily  on  aircraft  velocity,  which  does  not  vary 
significantly.  At  various  p,  the  closed  loop  phugoids  become 
unstable  with  a  natural  frequency  around  0.06 
radians/second.  This  instability  can  easily  be  compensated 
by  pilot.  This  is  also  true  for  rudder  failure  less  than  0 
degrees  as  seen  in  Figures  5-5  to  5-8. 

The  short  period  oscillations  (b)  in  the  open  loop  case 
presented  in  Figure  3-2  are  lightly  damped.  By  introducing 
the  controller,  their  damping  decrease  where  they  become 
unstable  for  values  of  p  between  -2.17  and  -1.28  degrees  as 
shown  on  Figures  5-3  and  5-4.  For  rudder  failure  of  -10 
degrees  the  instability  occurs  for  values  of  p  between  -2.66 
and  -0.7  degrees.  As  the  aircraft  moves  away  from  a  wings 
level  trim  condition,  a  decreases  and  the  bank  angle  becomes 
less  than  zero  degrees,  which  means  the  stability 
derivatives  related  to  or  are  affected.  Since  the  model  uses 
fixed  dynamic  stability  derivatives,  the  natural  frequency 
will  tend  to  decrease  and  the  damping  will  eventually 
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increase  (for  p  less  than  -2.17  degrees  for  a  0  degrees 
rudder  failure)  as  the  absolute  value  of  the  stability 
derivatives  related  to  a  decrease  [6:309].  For  a  rudder 
failure  of  less  than  -20  degrees,  the  short  period  roots  are 
pairing  with  the  controller  roots,  and  two  distinct  motions 
appear  as  shown  on  Figures  5-7  and  5-8.  The  sideslip  angle 
(8)  at  which  those  motions  occur  is  less  than  -7  degrees, 
which  is  outside  the  boundary  defined  earlier. 

Lateral  Motion 

In  Figure  4-3,  the  controller  has  lost  the  ability  to 
feedback  any  signal  to  the  rudder.  This  effect  is  seen  in 
the  dutch  roll  behavior  in  Figures  5-1  to  5-8.  Compared  to 
the  open  loop  case,  the  controller  basically  decreases  the 
dutch  roll  damping  since  it  has  lost  the  ability  to  feedback 
either  roll  rate  or  side  acceleration  to  the  rudder.  Since 
feeding  back  the  roll  rate,  in  this  case,  decreased  the  time 
constant  of  the  spiral,  this  also  had  the  effect  of 
decreasing  the  dutch  roll  damping.  As  8  moves  away  from  a 
zero  degrees  value,  the  real  part  of  the  dutch  roll  tends  to 
stay  stationary  while  the  magnitude  of  the  imaginary  part 
increases.  As  p  becomes  smaller,  the  dutch  roll  natural 
frequency  increases  and  its  damping  decreases. 

The  spiral  root  for  the  closed  loop  case  couples  with  a 
controller  root  to  create  a  pair  of  complex  conjugate  roots 
near  the  origin.  As  p  decreases,  the  roots  move  away  from 
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the  origin  and  create  a  lateral  short  period  oscillation 
before  returning  to  the  real  axis.  For  rudder  failures  less 
than  -10  degrees,  the  complex  conjugate  roots  break  on  the 
real  axis  where  they  become  real  as  seen  on  Figure  5-5. 

The  roll  also  couples  with  one  of  the  controller  roots. 
The  effect  encountered  from  that  coupling  is  only  apparent 
for  rudder  failure  less  than  -10  degrees  as  seen  on  Figure 
5-5.  A  complex  conjugate  root  is  formed,  moving  away  from 
the  origin,  which  produces  a  lateral  short  period 
oscillation.  As  /3  decreases,  the  roots  break  into  the  real 
axis  for  /3  less  than  -3.7  degrees. 

For  a  rudder  failure  less  than  -20  degrees,  the  real 
roots  located  on  the  real  axis  couple  to  form  a  short  period 
lateral  oscillation. 

For  case  A,  closed  loop  system  instability  occurs  for 
various  values  of  p .  This  instability  is  generated  by  the 
phugoid  root.  This  is  not  critical  since  the  natural 
frequency  of  the  phugoid  (around  .06  radians/second)  is  such 
that  the  pilot  can  easily  compensate  for  it. 

Discussion  of  Results 

Even  if  the  aircraft  can  be  trimmed  for  specific  rudder 
failure,  the  present  controller  modeled  in  Chapter  IV  puts 
restrictions  on  the  dynamic  stability  of  the  aircraft.  For 
the  open  loop  system,  the  eigenvalues  determined  for  both 
cases  gave  a  stable  system  for  p  less  than  -6  degrees.  Since 


5-4 


the  controller  loses  the  ability  to  use  the  rudder,  the 
control  laws  are  altered.  With  the  present  controller,  case 
A  is  stable  if  it  is  assumed  that  the  pilot  can  compensate 
for  the  phugoid  instability  generated.  For  case  B,  not 
taking  into  account  the  phugoid,  the  failure  permitted  with 
the  present  controller  is  as  followss 

Case  B,  0  degrees  rudder  failure:  the  aircraft  is 
stable  except  for  -2.2  <  /3  <  -1.3  degrees. 

Case  B,  -10  degrees  rudder  failure;  the  aircraft  is 
stable  except  for  -2.7  <  p  <  -0.9  degrees. 

Case  B,  -20  degrees  rudder  failure;  the  aircraft  is 
stable  for  -6.0  <  p  <  -3.3  degrees. 

Case  B,  -25  degrees  rudder  failure;  the  aircraft  is 
stable  for  -6.0  <  0  <  -4.9  degrees. 

In  Chapter  IV,  the  concept  of  controllability  was 
presented.  For  each  trim  condition,  the  open  loop  plant  and 
input  matrices  were  evaluated  for  controllability.  All  were 
controllable  with  the  available  input.  Therefore,  it  would 
be  possible  to  redesign  the  controller  in  order  to  prevent 
the  instability  of  the  closed  loop  system.  One  of  the 
important  factors  would  rest  in  the  control  power  stiil 
remaining  at  specific  trim  condition. 
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Summary 


For  both  cases,  the  phugoid  introduced  instability. 
Since  the  frequency  of  the  phugoid  is  low,  it  was  assumed 
that  the  pilot  would  be  able  to  compensate.  For  case  A,  the 
trim  regions  for  both  rudder  failure  types  were  dynamically 
stable.  For  case  B,  specific  trim  areas  within  each  rudder 
failure  types  are  to  be  avoided.  This  was  determined  using 
the  modified  controller  presented  in  figure  4-1  and  4-3. 
Reconfiguration  of  the  control  laws  could  be  done  to 
increase  the  trim  envelope.  The  only  limiting  factor  would 
be  the  control  power  still  remaining  after  the  aircraft  is 
trimmed . 

This  analysis  presented  an  appreciation  of  the 
different  types  of  motion  that  might  be  encountered  when  the 
aircraft  needs  to  be  trimmed  at  different  sideslip  angles 
due  to  a  rudder  failure. 
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LOCUS  OF  CLOSED  LOOP  EIGENVALUES 
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Figure  5-1  Case  A  0  Degrees  Rudder  Failure 
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LOCUS  OF  CLOSED  LOOP  EIGENVALUES 


Figure  5-2  Case  A  -10  Degrees  Rudder  Failure 
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REAL  AXIS 
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Figure  5-3  Case  B  0  Degrees  Rudder  Failure 
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Figure  5-4  Case  B  0  Degrees  Rudder  Failure 
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LOCUS  OF  CLOSED  LOOP  EIGENVALUES 


SIXV  AHVNIQVPMI 


Figure  5-5  Case  B  -10  Degrees  Rudder  Failure 
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(b)  SHORT  PERIOD  EIGENVALUES 


Figure  5-6  Case  B  -10  Degrees  Rudder  Failure 
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LOCUS  OF  CLOSED  LOOP  EIGENVALUES 


Figure  5-7  Case  B  -20  Degrees  Rudder  Failure 
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Figure  5-8  Case  B  -25  Degrees  Rudder  Failure 
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VI.  CONCLUSIONS  AND  RECOHHANDATIONS 


In  Chapter  I,  it  was  stated  that  this  research  would 
investigate  the  dynamic  stability  characteristics  of  an 
aircraft  which  has  sustained  failure  of  a  primary  control 
surface.  This  analysis  was  done  by  using  existing  data  to 
determine  specific  trim  conditions  and  to  evaluate  the  cross 
coupling  derivatives  to  be  included  in  a  linearized  aircraft 
model.  A  modified  controller  was  then  added  to  the  aircraft 
model  to  evaluate  the  closed  loop  response  of  the  system. 
Using  the  results  of  the  closed  loop  dynamic 

characteristics,  the  trim  regions  were  modified  to  represent 
the  damaged  aircraft.  The  following  paragraphs  provide  a 
summary  of  the  observations  and  conclusions  of  this 
research. 

Trim  Area 

Equi'i  jrium  analyses  was  performed  for  two  specific 
control  implementations.  The  first  used  the  actual  control 
surface  actuation  scheme  of  the  basic  F-16  aircraft.  The 
second  permitted  the  horizontal  tail  aileron  to  deploy 
independently  from  the  flaperon.  The  first  case  only  allowed 
trim  of  the  aircraft  for  a  rudder  failure  between  0  and  -10 
degrees.  The  second  case  allowed  trimming  the  aircraft  for 
the  worst  case  rudder  failure  analyzed  (-25  degrees). 


6-1 


Plant  Model 


When  the  equilibrium  areas  were  determined,  they  were 
used  as  input  to  the  linearized  equations  of  motion  of  the 
aircraft.  The  model  included  only  cross  coupled  static 
stability  derivatives.  The  aircraft  plant  matrices  were 
produced.  The  eigenvalues  of  the  plant  were  determined  to 
evaluate  the  dynamic  stability  of  the  system  for  each  trim 
condition.  It  was  observed  that  the  plant  matrices  were  all 
dynamically  stable  for  p  less  than  -6  degrees  at  the  given 
trim  conditions.  It  was  also  noted  that  for  all  trim 
conditions,  the  system  was  controllable.  Therefore,  it  is 
possible  to  improve  the  aircraft  response  even  if  it  has 
sustained  damage  to  its  rudder. 

Closed  Loop  System 

The  controller  developed  for  the  closed  loop  system 
took  into  consideration  the  damaged  rudder.  It  was 
interesting  to  note  that  the  controller  was  the  limiting 
factor  in  this  analysis  since  the  closed  loop  system  was 
dynamically  unstable  for  some  trim  regions  and  rudder 
failures.  The  results  also  showed  the  coupling  that  took 
place  when  trimmed  in  an  unsymmetrical  orientation.  As  p 
decreases,  the  roll  angle,  *,  decreases  which  introduces 
strong  coupling  between  the  lateral  and  longitudinal  motion. 
This  was  observed  by  analyzing  the  eigenvalues  of  the  closed 
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loop  system.  Satisfactory  aircraft  dynamic  response  would 
require  restructuring  the  control  laws  since  the  present 
controller  limits  the  aircraft. 

Recommendation 

During  the  course  of  this  thesis,  several  areas  of 
interest  emerged  which  would  provide  better  understanding  of 
the  problem  of  a  control  surface  failure.  They  are: 

1.  Tne  same  analysis  could  be  performed  without  using 
cross  coupling  stability  derivatives.  This  would 
provide  a  baseline  to  evaluate  their  effect  on  the 
aircraft  motion. 

2.  The  inclusion  of  the  cross  coupling  dynamic 
stability  derivatives  into  the  model  would  represent 
the  plant  more  accurately. 

3.  Development  of  a  new  controller  would  also  permit 
the  aircraft  to  sustain  damage  to  a  primary  control 
surface  by  using  greater  independence  of  each  control. 

4.  Only  rectilinear  flight  was  analyzed  in  this 
research.  Other  flight  conditions  could  also  be 
analyzed  to  determine  if  there  are  preferred 
trajectories  that  could  expand  the  trim  envelope. 

The  first  recommendation  could  easily  be  performed  by 
removing  the  cross  coupling  stability  derivatives  from  the 
actual  plant  model.  This  would  permit  the  evaluation  of  the 
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impact  of  the  different  derivatives  on  the  model. 

The  second  recommendation  would  require  the 
availability  cf  the  dynamic  cross  coupling  derivatives.  F^±. 
the  trim  conditions  evaluated  in  this  thesis,  it  would  be 
interesting  to  see  if  those  cross  coupling  derivatives  have 
a  noticeable  effect  on  the  observed  motion. 

The  third  recommendation  would  imply  the  redesign  of 
the  aircraft  control  laws  around  the  cases  investigated, 
which  is  a  major  task  since  several  trim  condition  were 
analyzed . 

In  the  last  recommendation,  combinations  of  aircraft 
trajectories  could  also  be  evaluated  together  in  order  to 
achieve  quas i - recti 1 inear  flight. 
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APPENDIX  A 


F-16  Layout,  Sign  Conventions,  and  Axis  Definitions 

Figure  A-l  shows  a  diagram  of  the  general  three-view 
layout  of  the  F-16.  Figure  A-2  defines  the  aircraft  axis 
systems,  and  the  angles  used  to  differentiate  between  them. 
Control  surface  deflection  conventions  are  also  shown  since 
definitions  for  positive  deflection  are  not  universal.  A 
graphical  representation  of  the  different  axis  system  is 
also  given. 


A-l 


-  WD 


*7  •  H  *  (4ft  * 


h A  •  *  #  <«fR  -  *rL> 


(1)  A  positive  control  fores  produces  e  oe$eclve  surfecs 
dafltctlon  and  eauaaa  a  poaltlva  aoaanc  about  aactl  aala. 

(2)  For  aactl  Individual  control  lurfaca,  tnlllnj-ad|a  down 
(laft)  La  poaltlva. 

O)  Laadln(-ad|a  flay  petition  la  oaaiurad  ttraatBrlaa.  *11 
octiar  poaltlona  art  oaaiurad  vich  ratpact  to  cha  bln«a* 
lino. 


Figure  A-2  Axis  System  and  Sign  Conventions 
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APPENDIX  B 


STATIC  AIRCRAFT  STABILITY  DERIVATIVES 


Included  in  this  appendix  are  the  functionnal  forms  of 
the  F-16  static  stability  derivatives  found  by  Zaiser  [5]. 
The  derivatives  are  used  in  the  trim  analysis  and  in  the 
state-space  model  of  the  aircraft.  Each  set  contain  the 
following  informations 

1.  The  control  surface,  i.e.  zero,  left  leading  edge 
flap,  etc. 

2.  The  force  or  moment. 

3.  The  correlation  between  the  experimental  data  and 
its  functional  form. 

4.  The  number  of  terms  in  the  polynomial. 

The  columns  of  the  data  file  contain  the  following 
information : 

1.  Number  of  the  polynomial  term 

2.  Power  on  the  a  term 

3.  Power  on  the  p  term 

4.  Power  on  the  6  term 

5.  The  coefficient  associated  with  that  term 
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APPENDIX  C 


Equilibrium  Area  Fortran  Code  Flow  Chart 

This  Appendix  contained  the  Flow  Chart  describing  the 
Fortran  code  used  to  perform  the  investigation  of  the 
equilibrium  area  for  case  A  and  B.  The  Flow  Chart  is  seen  on 
Figures  C-l  and  C-2. 
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C-l  Equilibrium  Area  Fortran  Code  Flow  Chart 
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Linearized  Equations  of  Motion 


Introduction 

In  order  to  gain  some  insight  into  the  nature  of  the 
dynamic  stability  characteristic  of  the  damaged  aircraft,  it 
is  necessary  to  derive  the  equations  that  governed  the 
aircraft  motion  about  some  nominal  condition.  This  Appendix 
will  take  the  general  equations  of  motion  that  are  developed 
by  McRuer  [6:203-232],  linearized  them  for  a  rectilinear 
flight  condition,  and  expands  the  force  and  moment  terms  to 
include  the  effect  of  coupling  that  are  generated  by  a  non 
zero  sideslip  angle  (p),  and  bank  angle  (*).  The  axis 
systems  used  in  this  appendix  are  the  aircraft  body  axis 
system,  and  the  stability  axis  system.  A  graphical 
representation  of  both  axis  system  is  shown  in  Figure  A-2 
of  Appendix  A. 

Equations  of  motion 

If  the  airframe  is  assumed  to  be  a  rigid  body,  the 
earth  to  be  fixed  in  space,  and  the  mass  and  mass 
distribution  of  the  aircraft  constant,  the  equations  of 
motion,  for  an  aircraft,  in  the  body  axis  system  are  given 
by  equation  D-l.  The  kinematic  equations  that  describe  the 
aircraft  attitude  are  given  in  equation  D-2.  All  force  terms 


D-l 


in  equation  D-l  incorporate  the  aerodynamic  and  the  thrust 
forces.  Assuming  that  the  aircraft  XZ  plane  is  a  plane  of 
symmetry  ( Iyy  and  Iyz  e<3ual  zero)  ,  equation  D-l  can  be 
simplified.  The  results  is  shown  in  equation  D-3. 


X  =  m  [  O+QW-RV  +  g  sin©  ] 

Y  =  m  [  V+RU-PW-g  cos  9  sin*  ] 

Z  =  m  [  ft  +  PV-QU-g  cos©  cos®  ] 

L  =  P  I  x  +  OR  (  Ij-Jy)  -  (  P  Q  +  R  )  IXZ 
+  (  P  R  -  0  )  Ixy  -  {  Q2  -  R2  )  IyZ 

m  =  o  Iy*  PR  (  ix  -  iz)  -  (  Q  R  +  p  )  Ixy 

+  (  P  Q  -  R  )  Iyz  -  (  R2  -  P2  )  Ixz 

N  -  R  I z  ■¥  P  Q  (  Iy  -  I x  )  -  (  P  R  *  0  )  Iyz 

+  (  Q  p  -  R  )  ixz  -  l  R2  -  02  )  ixy 


(D-l) 


*  =  P  +  Q  tan©  sin*  +  R  tan©  cos* 
6  =  Q  cos*  -  R  sin* 


♦  - 


cos* 

cos© 


sin* 

cos© 


(D-2) 


X  =  m  [  O+QW-RV+g  sin©  ] 

X  =  m  [  V+RU-PW-g  cos  0  sin*  ] 

Z  =  m  [  R+PV-ffl/-g  cos©  cos®  ] 

L  -  P  I x  +  Q  R  {  lx  -  Iy  )  -  (  P  0  +  R  )  Ixz 

H  =  0  Iy*  P  R  (  lx  -  Iz  )  -  (  P2  -  P2  )  Ixz 

N  .  R  IZ  +  P  0  (  Iy  -  Ix)  +  (  Q  R  -  P  )  Ixz 
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Linearized  Equations  of  Motion 

Since  equations  D-2  and  D-3  contain  products  of 
dependent  variables,  they  need  to  be  reduce  to  trackable 
form  where  the  total  motion  can  be  represented  by  an  average 
motion  (trim  condition  symbolized  by  the  subscript  0) 
representative  of  the  operating  condition,  and  a  dynamic 

motion  that  account  for  small  perturbations  about  the  mean 
motion  (symbolized  by  small  letter).  For  a  rectilinear 
flight,  the  average  motion  of  the  aircraft  about  it  center 
of  gravity  will  be  zero  as  well  as  all  acceleration  terms. 
Assuming  small  perturbations  about  the  trim  condition,  the 
linearized  equations  of  motion  can  be  written  in  the 
aircraft  body  axis  system  as 

dX  =  m  (  d+Wgq-Vgr+g  (cos 9q)  9  ) 

dY  =  m  (  V  +  £/0r-V0P-gr  cos00  cos*0* 

+  g  (  sin©0  sin*0)  9  ) 

dZ  =  m  (  (?  +  V0p-I/0g+g  (cos00  sin*0)  * 

+  g  ( sin00  cos*®)  9  )  (D-4) 

dL  =  p  1%  -  t  I  xz 

dM  =  q  IY 

dN  =  £  1%  -  P  I  XZ 
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and  the  linearized  kinematic  equations  as 

$  =  p  +  q  tan$0  sin#0  +  r  tan00  cos#0 
6  =  q  cos#0  -  r  sin*0 

(D-5) 

cos*0 
cos©0 

Assuming  that  the  thrust  vector  is  in  line  with  the  X  body 
axis  and  balance  all  other  trimmed  forces,  its  contribution 
to  the  perturbed  motion  can  be  neglected  [6:267].  Therefore, 
the  dX,  dY,  dZ,  dL,  dM,  dN  are  only  perturbed  aerodynamic 
forces  that  contribute  to  the  dynamic  of  the  aircraft.  Using 
small  angle  approximation,  it  is  possible  to  express  v  and  w 
in  term  of  /3  and  a  respectively.  The  resulting  equation  is 
shown  below  as 


f 

sin*0 

cos60 


(D-6) 


D-4 


Substituting  equation  D-6  into  D-4  and  rearranging  such 
that  the  time  derivative  terms  are  on  the  left  hand  side  of 
each  equations  will, lead  to 


a  = 


0  = 


r  11  -  Wq  g+  V0  r  -  g  {cos @0)  6 


—  -  Uq  r  +  V0  p 
m 


{  Uq2  +  Wq2  ) 


dr  = 


+  g  (cos©0  cos*0)  #  -  g  (sine®  sinig)  e  J 

{  n  r  ^ 

'  W 


-  Vq  p  +  Uq  q 


-  g  (COS00  sin#0  )  *  -  g  (sin00  cos*0)  e  J 


p  =  dL 


.  dM 

T-r 


jf  =  dL 


zz 

+  dN 

zxz 

zx  zz  -  zx/ 

\  ^ 

lx  ZZ  -  zx£ ; 

(D-7 ) 


Ixz 

+  dN 

Ix  ] 

/J <*Z-  ZXZ  t 

zx  zz  -  zxf 

i  =  p  +g  (  -ane0  sin*0)  +  r  {tan00  cos*0) 
d  =  q  cos *0  -  r  sin©0 


=  g 


sin#0) 

1  *  r 

f 

COS$0 

cose0  J 

COS00 

The  terms  that  need  to  be  determined  are  the  aerodynamic 
force  and  moment  that  are  included  in  equation  D-7. 

Aerodynamic  Forces  and  Moments 

The  representation  of  the  aerodynamic  forces  and 
moments  is  usually  made  in  the  stability  axis  system 


D-5 


(subscript  s).  The  aerodynamic  forces  and  moments  acting  on 
the  aircraft  in  the  stability  axis  system  are 

XS  =  -Drag(D)  =  -^LpV2SCD 
Yg  =  Sideforce(  Y)  =  -1  p  V2  S  Cy 
Zg  =  -Lift(L)  =  -  1  p  V2  S  CL 

(D-8) 

Lg  =  RollingMoment  =  -i  p  V2  S  Cj 
Mg  =  PitchingMoment  =  .1  p  V2  S  Cm 
Ng  =  YawingMoment  =  i  p  V2  S  Cn 


Since  the  equation  of  motions  are  written  in  the  body  axis 
system,  the  aerodynamic  forces  and  moments  need  to  be 
transformed  into  that  axis  system.  The  transformation  matrix 
that  puts  the  aerodynamic  forces  and  moments  into  the  body 
axis  system  is 


[BS] 


cosq  0  sine 
0  10 
-sina  0  cosa 


(D-9) 
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Transforming  the  aerodynamic  forces  and  moments  in  the  body 
axis  system  will  result  in 


X  =  X§  cos  a  -  Z  5  sincr 
=  j-i  p  V2  5  j  (  -  C  d  cos  a  +  CL  sina  ) 

=  1  p  V2  S  Cy 
Z  =  Xs  sina  +  Z §  cosa 

=  (  ^  p  ^  S  )  (  “  °D  Sinot  "  CL  cosa  ) 

L  =  Lq  cos  a  -  Nq  sina 

*  ( 4 » ^ s )  ( ci  cos  a  -  Cn  sina  J 
M  =  MS 

=  J  P  v2  s  cm 
N  =  Z,s  sina  +  N§cosa 

=  ( 4  p  ^ 5 )  ( ci  sina  +  Cn  cosa  J 


(D-10) 


Equation  D-10  can  be  expanded  to  determine  its 
dependence  with  the  perturbed  motion.  The  expansion  is  done 
using  a  Tailor  series  expansion  at  a  given  trim  condition. 


The  expansion  is  represented  by 


■L  -  ij)  + 

f  dF  ' 

MI 

k  y 

I0A1  + 

f  dF  ' 

MI 

k  > 

j  0X2  +  •  • •  + 

'  3F  ) 

M^J 

M  =  Mq  + 

r  dM  'i 

dAi 

k  / 

I0A1  + 

r  a/n 
d/^2 

k  ^  J 

1 0^2  +•  •  •  + 

'  a//) 
iMnJ 

0A./J 


0X/3 
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where  F  represent  the  forces,  M  the  moments,  A  the  dependent 
variables,  and  \  the  perturbation  from  the  trim  condition. 
The  dependent  variables  in  this  case  are 

A  =  u,  v,  w,  Q,  p,  q,  r,  6  (D-12  ) 


where 

u  =  ’/  cos/3  cosor 

v  -  V  sin/3  (D-13) 

w  =  V  cos/3  sina 

The  higher  order  term  have  been  eliminated  from  the 
expansion  to  keep  the  perturbed  equations  of  motion  linear. 
Taking  only  the  derivative  terms  in  equation  D-ll,  it  is 
possible  to  evaluate  the  dependence  of  each  aerodynamic 
forces  and  moments  in  the  body  axis  system  with  respect  to 
the  dependent  variables.  Since  the  aircraft  can  be  trimmed 
at  some  angle  a,  0,  and  #,  decoupling  of  the  linearized 
equations  of  motion  will  not  be  possible. 

Since  the  acceleration  in  the  X  body  direction  is  in 
general  negligible,  the  term  was  omitted.  Using  equation 
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D-ll  to  expand  equation  D-10  for  each  dependent  variables, 
this  will  give  the  contribution  of  the  forces  and  moments  to 
the  perturbed  motion.  Expanding  the  X  contribution  in  the 
body  axis  system  for  the  dependent  variable  u  will  give 


v  dD  n  3cos a 

=  Xn  =  -  — _  cosa  -  D  — _ 

u  du  du 


sin a  +  L 


3sin  a 


(D-14) 


where  the  derivative  of  each  components  of  equation  D-14  are 


dsing  =  dsina  fa  =  _  cosa 
du  da  du 


sina 


Uq  +  Wq‘ 


o 

dcosa  _  dcosa  da  _  sin  a 

-53-  *  ST-  53  =  ^ 2  +  „05 

iR  .  p  5  f  I  Cn  Rll  ♦  i  V2  w  HiE  ,  i  V2  ££? 

3u  K  2  "  3u  2  du  d  v  2  da  du 

+  £  ^2  dC£ 

2  3/3  au 

a  S  f  dCp 

=  -2__  2  Cp  cos/30  coscr0  +  Ve  — _  cos00  cosa0 

V0  l  a^ 

ac0  sina0  dCp  sin/30  cos/30  cosa0 


dCD  sina0 
l/0  _ —  —  -  *0 

Ju^+Wq2 


(Uq2  +  ^ 


_  A 

2  CL  COSP0  COSQ0  +  V0  — ^  COS00  C-OSO0 

au  1/0  a  v 


dCL  sina0 
0  3a  /  2 

_  i  2 

and  g  =  _  p  1/0 


3C£  sin/30  cos00  cosa0- 


(D-15 ) 
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Combining  all  the  terms  will  result  in 


*U  = 


g  s 
^0 


(  -2  Cq-  Cty  )  cos00  +  Vq 


sin^g  cos^0 


COS^O0 


Ju22  +  v 


^b2+*05 

(  (  ^  )  +  sin^0  cos^0  ) 

(-?  CL  -  c )  cos/30  j  sina0  cosa0 

^0 

+ - 

^02  +  V 


(  -  Ci^  -  Cp  )  sin 2a0 
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In  equation  D-16,  the  dimensional  derivative,  Xy  is 
expressed  in  terms  of  the  nondimensional  stability 
derivatives.  The  definition  of  each  nondimensional  stability 
derivatives  is  well  documented  in  McRuer  [ 6 : 292-293 ]. Also , 
the  addition  of  the  cross  coupling  stability  derivative 
terms  that  can  easily  be  related  to  the  both  longitudinal 
and  lateral  derivatives  for  definition.  The  only  difference 
with  McRuer  development  is  the  CiU  derivative  which  is 
define  as 


dCi 


~5v 


where  i  can  represent  lift,  drag,  side  force,  rolling 
moment,  pitching  moment,  and  yawing  moment. 

Similar  derivation  can  be  done  for  each  forces  with 
respect  to  each  dependent  variables  to  determine  the 
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expression  of  each  dimensional  derivatives  in  terms  of  the 
nondimensional  stability  derivatives  .  This  will  result  in 
the  following  equations 


YU 


q  s 


-(  -2  Cy  -  CYjj 


Vq  sin£0  cos00  ' 

- -  C  v. 

cosa0 

V 

W  ♦  *s 

- 

. *  .c,l 

sina0 
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Zu  = 


q  S 


(  -2  Cl  -  c iq  )  cospQ  +  Vq 


sin/30  cos^0 


2  +  Wo2 


■H 


COS 2O0 


/V  ♦ 


(  (  CL  -  )  -  sin^0  cosPq  Cq  ) 


-  {  -2  Cp  -  0^  )  cos/30  j  sincr0  cosa0 


^0 


/i/02  +  V 


(  CL  -  cn,  )  sin2a0 


2"':  „  5 


(D-18) 


XV  = 


q  s 
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^0  cosz/30 

{  -2  C p  -  CD  )  sin/30  -  _ - 1  C, 
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/i/02  +  Wtf 
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COSQ0 


sma0 
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q  S 
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g  g 
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1^0  COS  00 
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y0  COS ^00 
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Z  r  = 


q  S  b 
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\'CD, 
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X&i  =  q  5  |  -  C Di .  cosa0  +  CL sina0  j  (D-40) 

Yi>l  =  ?  5  Cirji  (D-41) 

^  =  g  5  /  -  C^.  sina0  -  .  cosa0  \  (D-42) 


Expanding  the  moments  contribution  with  respect  to  each 
dependent  variables  will  lead  to 
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Using  equations  D-16  to  D-69,  it  is  possible  to 
determine  the  contribution  of  the  total  perturbed  force  and 
moment  which  is 


d  X  =  Xu  u  +  X  v  v  +  X  w  w  +  Xp  V  +  Xff  Q 
+  Xp  p  +  Xg  q  +  Xr  r  +  X6.  6j 

d  Y  =  Yy  u  +  yy  v  +  w  +  yp  ?  +  Yft  & 

+  ypP  +  yg  g  +  yr  r  +  y6,  6^ 


(D-70) 
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It  is  now  possible  to  replace  equations  D-70  to  D-75 
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d  Z  =  Zj  u  +  Zv  v  +  Zww  +  Zp  V  +  Zft  # 

+  Zp  P  +  Zq  q  +  Zr  r  +  Z&i 

d  L  =  L(j  u  +  Lv  v  +  Ly/  w  +  Ly  V  +  Lft# 

+  L p  p  +  Lg  g  +  Lr  r  +  i6.  <5_j 

dM  =  My  u  +  MyV  +  Mff  w  +  Wp  y  +  0 

+  Wp  p  +  Mg  g  +  Wr  r  +  6  ^ 
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d  N  =  Ny  u  +  N  v  v  +  Nff  w  +  y  +  Nfl  ft 
+  N  p  p  +  N g  q  +  Nr  r  +  6  ^ 
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into  D-7  to  get  the  total  motion  of  the  aircraft  about  a 
trim  condition  which  can  be  represented  by 


X  =  A  x  +  B  u 


(D-76 ) 


where 


x  =  [  u  a  p  pgr#  0  +  ]T 

u  =  [&l  62  •  •  •  6n]T 

(D-77) 


and  6 n  represent  specific  control  surfaces  that  can  be 
controlled  by  the  pilot. 
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In  this  Appendix,  the  linearized  equations  of  motion 
for  an  aircraft  were  derived.  They  were  used  in  the 
computer  code  of  Appendix  E  to  get  the  F-16  plant  and 
control  variables  matrices  to  evaluate  the  dynamic  stability 


APPENDIX  E 


State  Space  Derivation  Fortran  Code 
This  appendix  describes  the  computer  code  used  to 
convert  the  stability  derivatives  for  an  aircraft  at  a 
specific  trim  condition  into  a  state  space  representation  of 
the  form 

X  =  A  X  +  B  U  (E-i) 


whe 

re  A  is 

the 

plant 

matrix , 

B  is  the 

control  matrix,  x  is 

the 

state , 

and 

u  the 

control 

.  Aircraft 

-specific  data  need  to 

be 

entered 

int 

o  the 

program 

as  well  as 

its  stability 

derivatives.  Cross-coupling  dynamic  stability  can  be  entered 
into  the  code.  Figure  E-l  is  a  symplified  flow  chart 
representation  of  the  code. 

If  more  than  one  trim  condition  is  input  into  the 
program,  the  output  for  the  A  matrices  is  a  file  where  each 
matrix  occupied  a  10x9  space.  Each  10x9  space  is  subdivided 
as  follows: 

a.  the  first  line  describe  the  aircraft  attitude  (0,  a, 

*,  9,  zeros) 

b.  the  remaining  9x9  space  represents  the  plant  matrix 

The  output  for  the  B  matrix  occupied  a  10x2  space.  The 

first  line  contains  £  and  or,  and  the  remaining  9x2  space  is 
the  control  variables  matrix  B. 


E-l 


E-2 


AFFEHDIZ  F 


CLOSED  LOOP  SYSTEM  EXAMPLE 


This  Appendix  contains  an  example  of  the  derivation  of 
the  closed  loop  system  eigenvalues  for  case  B. 

The  aircraft  states  are? 


X  = 


=  [  u  a  p  p 


u  = 


g  r  <j>  6  ♦  bFeff 

Jan a 


lr 


(F-l) 


The  open  loop  state  space  matrices  are  represented  by 


&  =  A  X  +  B  U 
X  =  C  X  +  D  U 


(F-2) 


where  A  is  the  plant  matrix,  B  the  control  input  matrix,  and 
C  and  D  are  the  measurement  matrices  associated  with  each  state  or 
control.  Matrices  A  and  B  were  developed  using  the  linearized  equations  of 
motion  of  Appendix  D  and  the  computer  code  developed  of  Appendix  E. 
Matrices  C  and  D  were  determined  using  equation  4-6  to  4-8c.  The  matrices 
are  presented  on  the  following  page. 
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A 


Colusns  1  through  6 


-1 . 31004-003 
-5 . 7800e-00S 
-2.  40004-006 
-  1 . 0000e-003 
1 • 93004-003 
3. 3100e-004 
0 
0 
0 
0 
0 


1 . 75004+001 
-8. 16004-001 
-1 . 47004-003 
2 . 95004+000 
-2 . 1600e+001 
-3. 1500e-001 
0 
0 
0 
0 
0 


1 . 09004+00 1 
2 . 9000e-001 
2 . 1 300e -00 1 
3. 70004+001 
7. 6500e+000 
1 . 4800e+001 
0 
0 
0 
0 
0 


-8 . 65004-004 
1 . 38004-002 
5 .45004-002 
-1. 44004+000 
-1.71004-003 
-1.41004-003 
1.00004+000 
0 
0 
0 
0 


-3.62004+001 
9 . 98004-001 
0 
0 

-5.51004-001 

0 

0 

1.00004+000 

-1.45004-002 


-9 . 23004+000 
0 

-9. 96004-001 
2. 18004-001 
0 

-4.73004-001 

0 

1 .45004-002 


0 


Columns  7  through  11 


-4 . 36004-005 
6 . 97004-004 
4.80004-002 
0 

-8.64004-005 

0 

0 

0 

0 

0 

0 


-3.21004+001 

-2.62004-003 

3.79004-005 

0 

3.25004-004 

0 

0 

0 

0 

0 

0 


0  -1.98004-001 
0  -3.57004-003 
0  -2.03004-003 
0  1.21004+000 
0  -2.89004-001 
0  5.17004-001 
0  0 
0  0 
0  0 
0  -2 . 0000e+001 
0  0 


-7 . 29004+000 
8.87004-004 
3.93004-004 
-4. 18004-002 
4.71004-002 
-4.42004-002 
0 
0 
0 
0 

-2.00004+001 


B 


D  = 


0  0 

0  0 

0  0 

0  0 

0  0 

0  0 

0  0 

0  0 

0  0 

20  0 

0  20 


0  0 
0  0 
0  0 
0  0 


C  = 


Columns  l  through  6 


0 

1.00004+000 

0 

0 

0 

0 

0 

0 

1 . 00004+000 

0 

0 

0 

0 

0 

1.00004+000 

1 . 97484-003 

6.74764+000 

-2 . 40954+000 

-2.72884-001 

-1.98934-001 

Columns  7  through  11 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1.37834-002 

5. 18084-002 

0 

-5.46224-002 

2.88384-003 

P-2 


<S>  ©  © 


The  state  space  system  that  represents  the  feedback 
paths  (Figure  4-4)  is  written  in  the  same  form  as  equation 


F-2  with  the  matrices  AK,  BK,  CK,  Dg.  The  original  system 
had  several  states  that  could  be  removed.  The  results  ares 

K  = 


-7.34974+000 
3 . 81614+000 
-7.6241e-016 
1. 3395e-016 


4 . 71516+000 
-3. 66994+000 
7.0566e-001 
-1. 15614-016 


-2.10754+000 
-1 . 0872e+000 
- 1 . 2343e+00 1 
2.690 le+000 


-2 . 9966e+000 
8.14334-001 
-5.1 9524+000 
1 .11214+000 


Bk  = 

0  6.085 14-003 
0  -8 . 6013e-003 
0  2.56164-002 
0  -4 . 4444e-003 


3 . 0575e-002  4 . 6053e-002 

-3 . 0824e-003  -3 . 1984e-002 
7.01224-002  -2 . 0804e -002 
-9.49724-003  5.0364e-003 


Ck 


-1.22704-014  1.06764-014  -1.99594-014  -8.33934+002 

0  0  0  0 


Dk  = 

-3.16804+000  0  0  -1.05894+000 

0  0  -1 . 2000e-001  0 

Assuming  that  no  inputs  are  feed  forward  ,  both  systems 
can  be  combined  together  using  equations  4-15  and  4-16  ,  to 
form  the  closed  loop  model  of  the  aircraft.  The  results  are 
shown  on  the  following  pages  using  the  representation  of 
equation  F-2  with  the  matrices  A^,  B^,  to 

represent  the  close  loop  system. 
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Columns  1  through  6 


-1. 3100e-003 

1 .  7500e+001 

1 . 0900e+00 1 

-8.65006-004 

-3 . 6200e+001 

-9 . 2300e+000 

-5 . 7800e-005 

-8 . 1600e-001 

2.  9000e-001 

1.38006-002 

9. 98006-001 

0 

-2  . 4000e-006 

-1 . 4700e-003 

-2. 13006-001 

5.45006-002 

0 

-9. 9600e-001 

-1.0000e-003 

2. 9500e+000 

-3. 7000e+001 

-1 . 4400e+000 

0 

2 . 1800e -00 1 

1 . 93006-003 

-2. 1600e+001 

7 . 6500e+000 

-1.71006-003 

-5.51006-001 

0 

3. 3100e-004 

-.3. 1500e-001 

1 . 4800e+001 

-1. 4100e-003 

0 

-4.73006-001 

0 

0 

0 

1.00006+000 

0 

0 

0 

0 

0 

0 

1 . 0000e+000 

1.45006-002 

0 

0 

0 

0 

-1.45006-002 

0 

4. 1824e-002 

2 . 06276+002 

-5. 1031e+001 

-5 . 77946+000 

-4 . 21316+000 

0 

0 

0 

0 

0 

2 . 40006+000 

0 

1 . 2017e-005 

7  .  1635e-002 

-1.46626-002 

4 . 4392e-002 

- 1 . 21056-003 

0 

-1.6986e-005 

-6 . 11206-002 

2 . 0725e-002 

-2.96366-002 

1 . 7110e-003 

0 

5.0587e-005 

2. 4297e-001 

-6. 1722e-002 

-2. 7794e-002 

-5. 0958e-003 

0 

-8. 7767e-006 

-3.94866-002 

1.07096-002 

6.24926-003 

8.84116-004 

0 

Columns  7  through  12 

-4. 36004-005 

-3.21006+001 

0 

-1. 9800e-001 

-7 . 29006+000 

0 

6 . 9700e-004 

-2. 6200e-003 

0 

-3.57006-003 

8. 8700e-004 

0 

4. 8000e-002 

3.79006-005 

0 

-2 . 0300e-003 

3.93006-004 

0 

0 

0 

0 

1. 2100e+000 

-«. 1800e-002 

0 

-8. 6400e-005 

3.25006-004 

0 

-2.89006-001 

4. 7100e-002 

0 

0 

0 

0 

5. 1700e-001 

-4.42006-002 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-2.91906-001 

1.09726+000 

0 

-2.11576+001 

6.10766-002 

2 . 4539e-013 

0 

0 

0 

0 

-2.00006+001 

0 

-8. 3869e-005 

3. 1526e-004 

0 

-3. 3238e-004 

1. 7548e-005 

-7 . 3497e+000 

1. 1855e-004 

-4. 4562e-004 

0 

4.69826-004 

-2.48046-005 

3.816 1 e+000 

-3. 5306e-004 

1.32716-003 

0 

-1.39926-003 

7.38726-005 

-7 . 6241e-016 

6.12556-005 

-2.30266-004 

0 

2.42766-004 

-1 . 2817e-005 

1.33956-016 

Columns  13  through  IS 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-2. 1352«-013 

3.99186-013 

1.66796+004 

0 

0 

0 

4.71 516+000 

-2 . 1075e+000 

-2.99666+000 

-3. 66996+000 

- 1 . 0872e+000 

8. 14336-001 

- 

7 . 05666-001 

-1.23436+001 

-5. 19526+000 

-1. 15616-016 

2.69016+000 

1 . 11216+000 

F-4 


0  0 

0  0 

0  0 

“i  0 

0  0 

0  0 

0  0 

0  0 

0  0 

20  0 

0  20 

0  0 

0  0 

0  0 

0  0 


Columns  1  through  6 


0 

1 ,0000e+00@ 

0 

0 

0 

0 

0 

0 

1 . 0000e+000 

0 

0 

0 

0 

0 

1 . 0000e+000 

1 . 9748e-003 

6.7476e+000  -2. 

.  4095e+000 

-2. 7288e-001 

-1.9893e-001 

Columns  7  through  12 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 . 3783s -002 

5. 1808e-002 

0 

-5. 4622e-002 

2.88386-003 

Columns  13  ' 

through  15 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0  0 
0  0 
0  0 
0  0 
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®  ®  ®  ®  <S>  ®  ® 


The  eigenvalues,  or  poles,  of  the  open  loop  plant 


are : 


0 

-6. 6025e-001+ 
-6. 6025e-001- 
-3.71 5 le -00 1+ 
-3.7 151e-001- 
-1 . 3977e+000 
-2 . 2807e-003+ 
-2 . 2807e-003- 
-2 . 8505e-002 
-2 . 0000e+001 
- 2 . 0000 e +00 1 


4. 7020e+000i 
4 . 7020e+000i 
4 . 0162e+000i 
4 . 0162e+000i 

6 . 4055e-002i 
6 . 4055e-002i 


The  eigenvalues  of  the  closed  loop  system  are  given  by 


0 

-2 . 0781e+001 
-  2 . 000 1 e+00 1 
-1 .  1034e+001  + 
-1 . 1034e+001- 
-8. 5100e-001+ 
-8. 5100e-001- 
-1 . 4483e-001+ 
-1 . 4483e-001- 
-9. 5788e-001+ 
-9. 5788e-001- 
-4. 2502e-003+ 
-4. 2502e-003- 
8. 5586e-003 
- 1 . 44  56e-001 


5. 9130e-001i 
5. 9130e-001i 
4. 6972e+000i 
4. 6972e+000i 
4. 1989e+000i 
4. 1989e+000i 
8. 1864e-001i 
8. 1864e-00 1 i 
5. 9495e-002i 
5. 9495e-002i 


And  the  transmission  zeros  of  the  close  loop  system  are 


-1 . 1250e+001- 
- 1 . 0000e+001  + 
-1 . 0000e+000 
-7. 5158e-008 
9. 1738e-017 


5. 1020e-016i 
1.0214e-016i 
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